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How ferromagnetic phases emerge in itinerant systems is an outstanding problem in quantum 
magnetism. Here we consider a repulsive two-component Fermi gas confined in a two dimensional 
isotropic harmonic potential and subject to a large Rashba spin-orbit (SO) coupling, whose single¬ 
particle dispersion can be tailored by adjusting the SO coupling strength. We show that the interplay 
among SO coupling, correlation effects and mean-field repulsion leads to a competition between 
ferromagnetic and non-magnetic phases. At intermediate interaction strengths, ferromagnetic phase 
emerges which can be well described by the mean-field Hartree-Fock theory; whereas at strong 
interaction strengths, a strongly correlated non-magnetic phase is favored due to the beyond-mean- 
field quantum correlation effects. Furthermore, the ferromagnetic phase of this system possesses a 
chiral current density induced by the Rashba spin-orbit coupling, whose experimental signature is 
investigated. 

PACS numbers: 03.75.Ss, 05.30.Fk, 75.70.Tj, 67.85.-d 


Introduction — Itinerant ferromagnetism represents 
an outstanding problem in many-body physics. Based 
on Stoner’s argument, a fermionic system in continuum 
may become ferromagnetic when the repulsive interac¬ 
tion strength exceeds a critical value [l|. In the context 
of a spin-1/2 Fermi gas, ferromagnetism means that the 
two spin species tend to phase separate to form spin do¬ 
mains as such a configuration obviously reduces interac¬ 
tion energy. Attempt to realize ferromagnetic state in 
repulsive Fermi gas was made by the MIT group in 2009 
0. Although some indirect evidences were present, spin 
domain formation was not observed. Later it was clar¬ 
ified that their system suffers from strong atom loss as 
the atoms tend to form tightly bound dimers, and ferro¬ 
magnetism was therefore not expected . From perhaps 
a more fundamental point of view, even if a repulsive 
Fermi gas is stable, it is not completely clear whether 
a ferromagnetic state will result. This is because the 
Stoner’s criterion is based on a mean-field argument, in 
which ferromagnetism arises once the mean-field repul¬ 
sion overcomes the kinetic energy. However, it has been 
conjectured that, under strong repulsive interaction, the 
Fermi gas may form a strongly correlated non-magnetic 
state [1, [HI that competes with the ferromagnetic state. 
Here the quantum correlation effects, neglected in the 
mean-field argument, play a more dominant role. There¬ 
fore whether itinerant ferromagnetic phases can exist in 
repulsive Fermi gases remains as an open question. 

In this Letter, we show that itinerant ferromagnetism 
can exist in a repulsive Fermi gas subject to spin-orbit 
(SO) coupling |6l4l2|. Itinerant ferromagnetism is a con¬ 
sequence of the interplay among kinetic energy, mean- 
field repulsion between the spin species, and quantum 
correlation effects. The key here is that the SO coupling 


significantly modifies the single-particle dispersion of the 
Fermi gas such that a ferromagnetic state emerges with¬ 
out the need of a very strong repulsive interaction. 


More specifically, we consider here a Rashba SO cou¬ 
pled repulsive two-component Fermi gas confined in a two 
dimensional (2D) isotropic harmonic potential. The trap¬ 
ping potential is necessary for any cold atom experiment 
as it provides atomic confinement. However, in the cur¬ 
rent situation, it plays an additional role: Together with 
the Rashba SO coupling, it produces a Landau level-like 
single-particle spectrum whose band flatness can be con¬ 
trolled by the SO coupling strength [iM3 , which as we 
will show is crucial for the existence of the itinerant fer¬ 
romagnetic phases in our system. If the system consists 
of a single spin-1/2 particle, the physics is well under¬ 
stood. The ground state is represented by a half vortex 
state [l^. An interesting feature of the system is that, 
under the limit of large SO coupling strength, the single¬ 
particle spectrum exhibits Landau level-like structure. In 
the case of an ensemble of spin-1/2 bosons, as previous 
works have shown 14, El , the near fiat band structure 
leads to intriguing spin textures and strongly correlated 
phases. Our current work tries to answer the question: 
What happens when we have an ensemble of repulsive 
spin-1/2 fermions? 


To address this problem and to elucidate the rela¬ 
tionship between ferromagnetism and interaction effects, 
we carried out a fully quantum mechnical exact diag- 
onalization (ED) calculation combined with a mean- 
field Hartree-Fock calculation. These two complemen¬ 
tary methods allow us to roughly divide the interaction 
effects into two parts: (1) it leads to a mean-field re¬ 
pulsion between the two spin species; (2) it builds up 
quantum correlations in the system. We will show that 
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FIG. 1. (Color online) Phase diagram in the g-\ plane. Here 
we consider 6 repulsively interacting spin-1/2 fermions con¬ 
fined in a 2D harmoinc trap, subject to Rashba SO cou¬ 
pling. g is the interaction strength, and is normalized to 
5 'mf = IM which is the mean-field critical interaction 
strength for a 2D Fermi gas without SO coupling. A is the 
dimensionless SO coupling strength. We can see that the 
phase diagram contains three phases: the weakly correlated 
non-magnetic phase, chiral magnetic phase and strongly cor¬ 
related non-magnetic phase. The dashed line represents the 
mean-field results which contains only two regimes: a non¬ 
magnetic phase below the dashed line and a ferromagnetic 
phase above the dashed line. 
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FIG. 2. (Color online) Left column: (al)-(a3) show the 
ground state angular momentum per particle J^/iV, entangle¬ 
ment entropy (EE) and ground state spin fluctuations (A^S^)^ 
as a function of interaction strength for A = 7, A = 6. The 
non-zero J^/iV indicates a magnetic ground state. Right col¬ 
umn: (bl)-(b3) show three representative single-particle oc¬ 
cupation number Um of state \m) for interaction strengthes 
marked by the yellow triangles in (al). 


the former favors ferromagnetism, whereas the latter has 
an opposite effect. The competition between them gives 
rise to the phase diagram shown in Fig. [U where a ferro¬ 
magnetic phase occupies a finite region in the parameter 
space spanned by the interaction strength and the SO 
coupling strength. 

The model — We consider a spin-1/2 Fermi gas, with 
atomic mass M and chemical potential confined in the 
x-y plane by an isotropic harmonic trap V{r) = 

(r = subject to a Rashba SO coupling Vsoc = 

^{Py<^x—PxO'y), where ax,y are Pauli matrices. The model 
Hamiltonian is given hy TL = TLo T 'Hint where 

no = J 

with ^ = (^>^,^^)^ being the atomic field operator, is 
the single-particle Hamiltonian, and 

Hint =9 j (2) 

with g > ^ decribes repulsive s-wave contact interaction. 
In what follows, we will adopt the trap units where the 
units for length and energy are given by Uho = ^/h/{MuJ) 
and huo^ respectively. Under this unit system, the in¬ 
teraction strength g has be units of huja‘^^ = h? /M. 
We also define a dimensionless SO coupling strength 
A = MAuho/^^- In the limit A ^ 1, the single-particle 
spectrum exhibits Landau level-like structure and the 


curvature of each Landau band is proportional to 1/A^, 
which provides a way to control the band flatness. Flat 
band structure will have two effects on an interacting 
many-body system: On the one hand, it may reduce the 
critical interaction strength for the ferromagnetic tran¬ 
sition according to the Stoner’s criterion. On the other 
hand, it makes quantum correlation more pronounced. 
Which of these two effects become more dominant deter¬ 
mines whether the system is ferromagnetic or not. 

ED results — The Landau level structure of the single¬ 
particle spectrum allows us to use the ED method to 
study a few-body system, where we restrict our calcula¬ 
tion to the lowest Landau level (LLL). The single-particle 
Hamiltonian Ho conserves the total angular momentum 
H, which is the sum of the orbital and the spin angu¬ 
lar momentum. Single-particle states in the LLL can be 
labeled by a single quantum number |m), whose total 
angular momentum is H = m + 1/2, and whose energy 
(apart from a constant) is approximately m{m + 1)/A^ 
|l3l4l5|. A set of such states form a Fock space basis, 
upon which the total Hamiltonian can be expanded 0- 
For details, see Supplemental Material [T^ . 

We present our ED results for a system of A" = 6 
fermions. Under the total Hamiltonian H, H of the whole 
system remains as a good quantum number. Figure [21(a) 
displays H of the ground state as a function of the inter¬ 
action strength g for a fixed SO coupling strength A = 7. 
As one can see, Jz = 0 for small g^ becomes finite for 


2M 


— /i + Vsoc + U(r ) 
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FIG. 3. (Color online) (a) shows the spin texture of the ferro¬ 
magnetic state and (b) shows the total number density (back¬ 
ground color) and the chiral current density (arrows) of the 
ferromagnetic state. 


intermediate and vanishes again at large g. The single¬ 
particle occupation number rim = {cilri^rn)^ with the 
annihilation operator associated with state |m), for the 
three representative cases are plotted in Fig. [2jbl)-(b3) 
as a function of m. For a weakly interacting system, as 
shown in Fig. [2](bl), interaction induces a few particle- 
hole excitations near the “Fermi surface”. However, the 
ground state still preserves the time reversal symmetry, 
i.e., Tirn = n-m-i‘ The density profiles for the two spin 
species are identical. One can also calculate the local spin 
vector s{r) = (T’t'(r)aT(r)) and show that it vanishes ev¬ 
erywhere. Hence the state is a non-magnetic state. 

At intermediate interaction strength, as shown in 
Fig. [2](b2), the ground state breaks time reversal sym¬ 
metry with rim 7 ^ the two spin species possess 

non-overlapping density profiles, and non-vanishing local 
spin vector s{f) emerges, see Fig. [31(a). This indicates 
that the state is a ferromagnetic state. Furthermore, we 
calculated the current density of this state. With Rashba 
SO coupling, the current density is given by 

jir) = + Xzxs, (3) 

m 

where - 4>\ry4>m]/‘2: comes from the 

orbital motion where c^m represents the wave function of 
the single-particle state \m). Due to the the time reversal 
symmetry of the Hamiltonian, we have 
In the non-magnetic state, rim = and s = 0, both 

terms on the right hand side of Eq. vanish. However, 
for the ferromagnetic state, they are both finite and leads 
to a chiral current as shown in Fig.|3l(b). As a result, we 
call the magnetic state chiral ferromagnetic. 

At large interaction strength, as shown in Fig. |2](b3), 
the time reversal symmetry is restored, and once again we 
have s{f) = 0 and j{f) = 0 as in the weakly interacting 
regime. The fluctuations of rim indicates that this non¬ 
magnetic state is strongly correlated. To quantify the 
quantum correlation and fluctuation, we calculated the 
entanglement entropy (EE) of the system [l^, and the 
total spin fluctuation {ASz)‘^ = {S‘^) — {Sz)‘^^ and plot 


them as functions of g in Eig.[2](a2) and (a3), respectively. 
Both EE and {ASz)^ for the large interaction regime are 
significantly higher than those in the other two regimes. 


With the above results and similar calculations for 
other SO coupling strengths, we can present the phase 
diagram as shown in Eig. [H Eor A < 13, there exists 
a window of ferromagnetic phase at intermediate values 
of As ^ increases from zero to a lower critical value 
(represented by the red solid line with filled circles), the 
weakly correlated non-magnetic state becomes ferromag¬ 
netic. Note that this lower critical value is much smaller 
than = 27r/i^/M, the mean-field ferromagnetic criti¬ 
cal interaction strength of a 2D homogeneous Eermi gas 
without SO coupling [illii. This can be understood 
from the Stoner’s argument and the fiat band single¬ 
particle spectrum. In Ref. 0 , it was shown that the 
critical interaction strength for ferromagnetic transition 
in a repulsive Eermi gas can also be reduced by adding 
a weak optical lattice, as the lattice potential helps to 
quench the kinetic energy. The essential physics here is 
similar to our situation. Howver, as g further increases 
to an upper critical value (represented by the blue solid 
line with empty squares), the ferromagnetic state gives 
its way to a strongly correlated non-magnetic state. As 
A increases, i.e., the single-particle band becomes flatter, 
this window of ferromagnetic phase shrinks quickly, and 
eventually vanishes for A > 13. At such large SO cou¬ 
pling strength, the single-particle band becomes so fiat 
that a very small interaction strength gives rise to strong 
correlations that disfavor the ferromagnetic state. 


To examine the finite-size effect, we made ED calcula¬ 
tions for A" = 4, 6 , 8 , and found that with increasing A, 
the quantum correlation effect is somewhat weakened. 
Eor fixed A, the lower critical interaction strength at 
which the weakly correlated non-magnetic phase changes 
to the ferromagnetic phase is not very sensitive to A, 
while the upper critical interaction strength at which 
the ferromagnetic phase becomes the strongly correlated 
non-magnetic phase increases with A. Eurthermore, the 
critical SO coupling stength at which the ferromagnetic 
window vanishes also increases with A. As a result, the 
ferromagnetic regime in the g-X phase diagram is en¬ 
larged with increasing A. 


Hartree-Fock results — If it is the correlation effects 
that destroy the ferromagnetic state, then one should 
not expect this to occur in a mean-field theory, which 
neglects quantum correlation. To examine this, we now 
turn to a mean-field Hartree-Eock (HE) calculation. Un¬ 
der the HE theory, the many-body wave function takes 
the form: = ^ 7 ^ Ep(-l)^ 0 i(n)</> 2 (f^ 2 )--</>Jv(nv), 

where P represents permutations, and are single- 
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particle orbitals that satisfy the following HF equations: 


+ iX {—dyCFx + dxCFy) + 

+ fpa{^=^a 0a, (4) 


where n(r) = X:a=i l0aP and m (r) = Y!1=i 0l^0a are 
local density and spin vector, respectively. 

We numerically solve the HF equations self- 
consistently without invoking the LLL approximation (for 
details, see Supplemental Material [l^ ). In Fig. lUwe plot 
density profiles from this calculation. Here we also take 
N = 6 and A = 7 in order to make comparisons with the 
ED results. However, we also performed HF calculations 
up to = 200 and found no qualitative differences from 
the N = 6 results presented here. For small interaction 
strength g = O.OlS^^p [Fig. HJ^a)], both HF and ED tell 
us that the state is non-magnetic with identical density 
profiles for both spin species. Eurthermore, the results 
from the two theories agree with each other very well. 
At ^ = O.OS^f^p [Eig. |4l[b)], ED predicts a non-magnetic 
state, whereas HE indicates that the system already en¬ 
ters the ferromagnetic regime. In fact, HE calculation 
predicts a critical interaction strength ^hf ~ 0.027^mf 5 
while the corresponding critical interaction strength for 
ED is gEB ~ 0.05^mf- That ^ed > ^mf can be at¬ 
tributed to the fact that the quantum correlation in the 
ED calculation disfavors the ferromagnetic phase. At 
g = 0.173^^p [Eig. [4](c)], HE and ED agree with each 
other again, both predicting a ferromagnetic state. At a 
large interaction strength g = 0.2345f^p [Eig. [IKd)], dis¬ 
crepancies arise between the two calculations again: ED 
predicts a non-magnetic state, while HE gives a ferromag¬ 
netic state. In fact, as we have expected, for g > guF^ 
HE always predicts a ferromagnetic state. In contrast, 
our ED calculation shows that for sufficiently large g, 
strong correlation destroys the ferromagnetic state. Eor 
the parameters used in Eig. 01 ED shows that ferromag¬ 
netic phase only exists for 0.05 < g/g^E ~ 

In the phase diagram of Eig. [H the dashed line rep¬ 
resents ^hF: which separates the phase space into non¬ 
magnetic (below the dashed line) and ferromagnetic 
regimes (above the dashed line), ^hf decreases quickly 
as A increases (which can again be understood as due to 
the band flattening), but never terminates as in the case 
of ED. To further demonstrate the effects of quantum 
correlation, we plot in Eig. [5] the energy as a function of 
interaction strength at A = 7. Eigure[5](a) shows how the 
total energy the kinetic energy Ekin and the inter¬ 
action energy E[nt from the ED calculation change as g. 
As g increases. Eg keeps increasing monotonically, while 
Eint decreases at both phase transition points (shown 
by the vertical lines) at the cost of increasing Ekin- In 
Eig. [51(b) we plot the ratio of the interaction energy, which 
is simply g times the density-density correlation between 





FIG. 4. (Color online) Density profiles of each spin species for 
different interaction strengths with A = 6 and A = 7, by both 
ED (red solid lines: thick lines for spin up and thin lines for 
spin down) and HF (black dashed lines: thick lines for spin 
up and thin lines for spin down) methods. 


the two spin species integrated over all space, from the 
ED and the HE calculation. As it shows, in the weakly 
correlated non-magnetic and the ferromagnetic regimes, 
the ED and the HE results are comparable to each other. 
By contrast, in the strongly correlated regime, the in¬ 
teraction energy from the ED calculation is signaficantly 
lower than that from the HE calculation. This clearly 
shows how the system develops nontrivial quantum cor¬ 
relations such that, even though the density profiles of 
the two spin species completely overlap with each other, 
the joint probability of finding two unlike spins at the 
same position is strongly suppressed. 




FIG. 5. (Golor online) (a) ED results for the ground state en¬ 
ergy Eg, the kinetic energy Ekin, and the interaction energy 
Eint as functions of interaction strength, (b) The ratio of the 
interaction energy from the ED calculation and that from the 
HF calculation. The two vertical lines separate the parame¬ 
ter space into three phases according to the ED calculation: 
from left to right, we have the weakly correlated non-magnetic 
phase, the ferromagnetic phase, and the strongly correlated 
non-magnetic phase. Here A = 6 and A = 7. 

Einally, we propose an experimental procedure to de¬ 
tect the chrial current associated with the ferromagnetic 
state. The procedure goes as follows: Eirst the ground 
state (either magnetic or not) is prepared. Then the 
harmonic trap is suddently distorted from isotropic to 
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results for the repulsive Fermi gas remain essentially the 
same under Dresselhaus coupling 0- 
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FIG. 6. (Color online) Time evolution of the atomic cloud 
after a suddent trap deformation. At t = 0, the trapping fre¬ 
quency along the y-axis is suddently changed from uj to 3.16cc;, 
while that along the a:-axis remains at uj. The upper (lower) 
panel shows the dynamics of a non-magnetic (ferromagnetic) 
state. Here A = 6, A = 7 and tq = 1/uj. 


anisotropic. For a non-magnetic state, as shown in the 
upper panel of Fig. [7| obtained from a time-dependent HF 
calculation 0 , this induces a quadrupole mode. By con¬ 
trast, for an initial chiral ferromagnetic state, the whole 
cloud also undergoes an angular rotation, analogous to 
the scissors mode in a condensate with vortices [20|. 

Conclusion — In summary, we have shown that how 
the Landau level-like band structure of a 2D Rashba SO 
coupled Fermi gas, with a controllable band flatness, can 
be exploited to exhibit itinerant ferromagnetism. The 
near-flat band structure dramatically reduces the critical 
interaction strength required for the ferromagentic phase 
transition. We employed two complementary methods, 
the fully quantum ED method and the mean-field HF 
method, to investigate this problem. Our calculation 
elucidates the interplay between the mean-field repulsion 
and the quantum correlation effects, and shows that the 
former favors while the latter tends to destroy ferromag¬ 
netism. The emergence and disappearance of the fer¬ 
romagnetic phase result from the competition between 
these two factors. We have also shown that the ferro¬ 
magnetic phase in our system is accompanied by a chiral 
density current resulting from the SO coupling. This chi¬ 
ral density current and the spin texture that characterize 
the itinerant ferromagnetic state can be readily detected 
using today’s cold atom techniques. We hope our work 
may open new avenues of research in both SO coupling 
and itinerant magnetism in cold atoms. 


Finally, we comment that Dresselhaus SO coding has 
recently been realized by the Shanxi group [2l|. The 
single-particle spectrum of a harmocally trapped 2D spin- 
1/2 particle remains exactly the same if the Rashba SO 
coupling is changed to the Dresselhaus coupling. Our 
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Supplementary Materials for “Itinerant chiral ferromagnetism in a trapped Rashba spin-orbit coupled Fermi 

gas” 


In this Supplementary Materials, we provide more technical details of the main manuscript. 


Exact diagonalization scheme 


For a single spin-1/2 particle with Rashba spin-orbit (SO) coupling confined in a two dimensional (2D) isotropic 
harmonic trap, the single-poarticle Hamiltonian is given by Ho in the main text. The spatial wave function of the 
eigenstates take the form: 


(^ ) 


fnm ('^) 

9nm(.‘r)e^* 


^imcf) 


^n,—m—1 (^ ) 


9nm{r)e 




n = 0,1,2, 


m = 0, dzl, dz2, ... 


(5) 


which form a degenerate time reversed pair with eigenenergies €n,m = In the limit that the dimensionless 

SO coupling strength A ^ 1, the eigenenergies (apart from a constant) take the following approximate form: 


^n,m — ^n,—m—1 


n + 


m{m -h 1) 
A2 


huj . 


For a few-body system with weak interaction and small particle number, the Hilbert space is limited to the lowest 
Landau level (LLL) which is specified by the quantum number n = 0. We introduce a cutoff m* which further reduces 
the Hilbert space to that with —m* — 1 < m < m*. The value of m* is determined by specific values of N and g. Given 
N fermionic particles filled to M = 2m* + 2 single particle states, we obtain totally Fock states. Due to the 

rotational symmetry of this system, we are able to divide the full truncated Hilbert space into several independent 
subspaces with fixed total angular momentum which considerably reduces the dimension of the 

Hamiltonian that needs to be diagonalized. 

Next we present the main steps for the ED scheme for specific subspace with and particle number N. The Fock 
states are denoted by \pi) = = 1,2, ...,D with the convention mi < m 2 < ... < mjsf- For later use, 

we associate each occupied single particle-state rrii with a number (for example, = 1). The single-particle 
part of the Hamiltonian Ho is diagonal under this basis: 

N 

{Pi\'^o\Pj) = . (6) 

a=l 

where = 6 ^= 0 ,m- Under the same basis, the diagonal matrix elements of the interacting Hamiltonian Hint take the 

form: 

(PilHintbi) = I j d'^rip^ - , (7) 

where p{r) = {pi\p\pi) and s = {pi\^Pi) represent the local density and spin vector, respectively. This diagonal matrix 
elements can be regarded as the mean-field Hartree-Fock interaction energy associated with the Fock state \pi). The 
non-diagonal matrix elements of Hint are non-vanishing only between two Fock states that differ by two single-particle 
states, say \p) = ...aj^...al ...| 0 ) and \q) = ...a^...aj...| 0 ) with the constraint m + n = /c + /: 

{q\Uint\p) = J , 

where cr | denotes the wave function of the single particle state in the lowest Landau level. The non-diagonal 
part of Hint builds up correlations between different Fock states. It mixes Fock states with different spin polarization, 
therefore tends to suppress the magnetic phase. 
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Calculation of entanglement entropy (EE) 


Entanglement measure is useful to analyze correlation properties of the ground state. We calculate EE in the 
following way. We first divide the system into two subsystems (denoted as A and B) and then analyzing the reduced 
density matrix in one of the subsystems. In our system, the subsystems can be distinguished by the single particle 
angular momentum jz = m + the Subsystem A includes all the positive jz states, while the Subsystem B includes 
all the negative jz states. The total ground state density matrix is given by p = \G){G\ with \G) denoting the ground 
state. By the standard procedure, we trace out the Subsystem B to find the reduced density matrix for Subsystem A: 

JC __ 

where jc = mc^l/2 denotes a finite-size cutoff of this system. The eigenvalues of the reduced density matrix pf give 
rise to the entanglement spectrum = — Inp^. Eor pure Eock state without any correlation, there will be only one 
non-zero eigenvalue pf = l and all the others equal to zero. Therefore, we can observe only one point with ^ 0 
and other points 1 in the entanglement spectrum for less correlated ground state. While for strongly correlated 
ground state, the entanglement spectrum has a broad and fiat structure. We can further calculate the ground state 
EE by EE = —trp^ Inp^ = — pf Inpf. We will find EE ^ 0 for less correlated ground state while EE > 1 for 
strongly correlated ground state. 


Hartree-Fock equation for trapped spin-orbit coupled Fermi gas 


For weakly correlated states, the mean-field Hartree-Fock (HE) approximation captures the key physics. The HE 
approximation neglects quantum correlations of the state by assuming: 


Thf 




(9) 


where P denotes all permutations, and are orthonormal single-particle orbitals to be determined. With this 
assumption, we can obtain the HE Hamiltonian as follows (adopting the trap units): 


i^HF 


J 


■iv^ + iX {-dya^ + d^ay) + + |n(r) - |m(r) 




( 10 ) 


where the constant terms (p/4) j dr[m{rY — n(r)^] has been dropped, and n(r), m(r) are respectively the averaged 
local density and spin vecotr: 


N N 

n{r) = , m(r) = ^4(r)CT • (11) 

q;= 1 a=l 

To manipulate the interaction term in a spin rotational invariant way, we have rewritten the interaction term as: 
§ f dr (n^ — ^) in the above calculation. The single-particle wave functions (ja with a = 1,2, ...,A/' satisfy the HE 
equations: 


+ iX{-dya^ + d^dy) + E + |n(r) - |m(r) 




( 12 ) 


which, together with Eq. (HID, form a closed set and can be solved self-consistently. 

In our calculation, 0 q,’s are expanded onto the single-particle eigenstates defined in Eq. ([5]): = 

Snmfhat in our HE calculation, we do not restrict to the LLL. So we have to introduce a 
cutoff Nc for quantum number n, in addition to the cutoff for quantum number m. Under this expansion, the HE 
equations take the form: 


E(< 

n2 = l 


nim'-'nin2 


I A ^ n. 



'^a-,n2m 




(13) 
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where (r), = J • m(r). Due to the rotational symmetry, m is a 

conserved quantum number. The Hartree-Fock wave function Thf would be obtained through iteratively solving the 
above equations until self consistency is reached. While the ED calculation can only hand a few particle numbers (up 
to 8 in our calculation), we have done HE calculations up to 200 particle number. Erom our calculation, we found that 
the mean-field critical interaction strength at which the non-magnetic state changes to ferromagnetic state roughly 
scales as ^hf cx: N /\^. 


Time-dependent Hartree-Fock theory 

To study the dynamics, we extend the HF calculation to time-dependent situation. The time-dependent Hartree- 
Fock equations take the form: 

+ *A {-dyax + dx(Jy) + V{f,t) + J [n (r, t) - m (r, t) ■ (!>„ (r, t) = idt(l>a {r, t) , (14) 

where we have assumed that the trapping potential V{f,t) is time-dependent. The initial wave function at t = 0 is 
taken as the ground state wave function under E(r, t = 0). The orthonormality of the single-particle orbitals (j)a{r,t) 
at time t is guaranteed by the unitary time evolution. The local density and spin vector n and m are still given by 
Eqs. ([TT]) with the explicit time dependence. 
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FIG. 7. (Color online) (Color online) Time evolution of the atomic cloud after a suddent trap deformation. At t = 0, the 
trapping frequency along the y-axis is suddently changed from uj to V^<^, while that along the a:-axis remains at u. The upper 
(lower) panel shows the dynamics of a non-magnetic (ferromagnetic) state. Here to = 1/ct;, A" = 200 and A = 20. 

In the main text, we presented results for 6 fermions. In Fig. [71 we show the similar time evolution for 200 fermions. 
The key physics remains the same: Under a sudden deformation of the trapping potential, the non-magnetic state 
exhibits quandrupole oscillation, while the chiral ferromagnetic state exhibits an additional rotation analogous to the 
scissors mode. 


Mapping to Dresselhaus SO coupling 

In our calculation, we have taken the SO coupling to be of Rashba form: Vsoc = ^{PyO'x — PxO'y). The results can 
be easily generalized if the SO coupling is of Dresselhaus form: V|q = ^{PyCfx PxO'y). The system with Dresselhaus 
SO coupling can be mapped to a system with Rashba SO coupling through a unitary transformation in spin space: 
U = icFx^ under which the Pauli matrices are transformed as 

<^x ^ ^ ^y ^ ~^y 5 —(^z 5 

and the Rashba SO coupling is then transformed to the Dresselhaus form. The s-wave interaction is spin SU{2) 
invariant and will not be changed under the above unitary transformation. So, all the results achieved in our main 
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text hold for Dresselhaus SO coupling case after this spin space transformation. For example, the single-particle wave 
function is obtained by the single-particle and many body energy spectra are unchanged. The main 

difference is the change of the ground state spin texture, where the local spin transforms as: 


as schematically shown in Fig. [HI The density current operator coming from the Dresselhaus SO coupling is also 
modified to be = A(sy, 5^,, 0), which keeps the density current in the ground state invariant (Fig. [8]). 



(b) 



FIG. 8. (Color online) The schematic plot of the spin polarization (black arrow) and density current (red arrow) for (a) Rashba 
SO coupled system and (b) Dresselhaus SO coupled system, respectively. 







